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Nanoparticles with multiple ligands have been proposed for use in nanomedicine. The multiple 
targeting ligands on each nanoparticle can bind to several locations on a cell surface facilitating both 
drug targeting and uptake. Experiments show that the distribution of conjugated ligands is unex- 
pectedly broad, and the desorption rate appears to depends exponentially upon the mean number 
of attached ligands. These two findings are explained with a model in which ligands conjugate to 
the nanoparticle with a positive cooperativity of « 4fcT, and that nanoparticles bound to a surface 
by multiple bonds are permanently affixed. This drives new analysis of the data, which confirms 
that there is only one time constant for desorption, that of a nanoparticle bound to the surface by 
a single bond. 



A dendrimer is a branched polymeric nanoparticle with 
the topology of a Cayley tree [1]; see Figure[l] We will be 
concerned here with dendrimers with a radius of about 5 
nm having w 100 termini that can be functionalized by 
the conjugation of various endgroups and ligands. These 
terminal groups can be varied to tune solubility proper- 
ties, and different ligands can be used to target and treat 
various cell pathologies. [21 E] [4]. 

Targeting ligands can be used to enhance the binding 
of the nanoparticle to specific receptors [5]. For example, 
epithieal cancer cells are known to overexpress folic acid 
receptors, so that folic acid attached to the dendrimer 
should target epithelial cancer, allowing chemotherapy 
agents also attached to the dendrimer to have high speci- 
ficity [6j. For this application, it is important to under- 
stand the statistical distribution of the number of at- 
tached ligands and how this distribution affects binding 
to the cell surface. That is the subject of this paper. 
As we will see, this distribution is broad so that the 
fluctuation of the number of ligands is comparable to 
its mean. These large fluctuations are characteristic of 
physics at the nanoscale; this effect is often neglected. 
Further, the chemical reactions in question are always 
far from equilibrium. Analysis of the nanoparticle prod- 
uct from these reactions is in this case unique because 
detailed data is available on the ligand number distri- 
butions. This data allows us to explore effects, such as 
cooperativity in ligand conjugation during synthesis and 
multivalent enhancement of binding affinity, that would 
otherwise be unaccessible. 

In our experiments measuring the distribution of lig- 




FIG. 1: (Color Online) A dendrimer bound to a surface. 
The dark spheres represent folic acid molecules and the light 
spheres folate binding protein, the folic acid receptor, on the 
surface. Molecular dynamics simulation by C. Kelly. 



and attachment [7j we conjugated varying amounts of 
the ligand 3-(4-(prop-2-ynyloxy)phenyl)propanoic acid to 
the surface primary amines of a poly(amidoamine) den- 
drimer (G5 PAMAM; (NH2)iio). This ligand was chosen 
because its binding properties and steric constraints are 
similar to folic acid, and because it is amenable to separa- 
tion by High Pressure Liquid Chromatography (HPLC). 
Ligation takes place in a solution with n ligands available 
per dendrimer; the conjugation is by random attachment, 
and we assume that all the ligands attach. Let C„ be the 
distribution of ligands on the dendrimers. C„ is measured 
by HPLC. We found that C„ is very broad, in fact broader 
than a Poisson distribution; see Figure [2j We attribute 
this effect to cooperativity, i.e., binding one ligand makes 
it easier to bind more ligands. 
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FIG. 2: Distributions of dendrimers with n bound ligands, for 
average ligand number n of (a) 0.9, (b) 3.7, (c) 5.8, and (d) 
13.9. The data points are HPLC data from ligation exper- 
iments that are chemically similar to the folic acid ligation. 
The dashed line is the Poisson distribution and the solid line is 
the result of Monte-Carlo simulations with E = — 3.7±0.1kT. 



The dendrimer was experimentally determined to have 
approximately 110 free sites at the end of its (roughly) 
spherically arranged branches. This is modeled as a 
11x11 triangular lattice with periodic boundary condi- 
tions. We use a kinetic model of cooperative ligation 
with two parallel attachment paths. In the free attach- 
ment path, a ligand attaches with a free energy barrier 
Eq = 0. In the cooperative path, a ligand attaches with 
free energy barrier i? to a site which has ligated neigh- 
boring sites. 

For large n, two factors must be considered. First, 
when a reaction site is proximate to a previously ligated 
site there is a possibility for a catalytic enhancement of 
the reaction rate. The reaction occurs at the amide group 
of the ligand-nanoparticle bond, so the presence of mul- 
tiple neighbors can increase the probability of a ligand 
attaching, say by orienting it properly [S]. Second, mul- 
tiple ligand neighbors crowd the site, sterically hindering 
ligation. The data suggest that the hindrance is so strong 
that ligation does not occur when there is more than one 
ligand neighbor. We write for the rate of attachment 
at a target with with n; ligated neighbors and n/ free 
neighbors, where ni + Uf =6: 



R = LOoL(t)S{ni){nit 



-E/kT 



-nf)/{nf + ni) 



(1) 



Here is a molecular time scale, L(t) is the free ligand 
concentration, and S{ni) is a steric hindrance term that 
is equal to one if n; = {0, 1} and is zero otherwise. 
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FIG. 3: Distributions of dendrimers with 68% of sites acety- 
lated and with n bound ligands, for average ligand number n 
of (a) 1.0, (b) 1.9, (c) 3.4, and (d) 16.0. The data points are 
HPLC data. The dashed line is the Poisson distribution and 
the solid line is the result of best-fit Monte-Carlo simulations 
with E = -5.0 ± 0.2kT. 



We implemented a continuous-time, rejection-free 
Monte-Carlo simulation of the model with N = 10000 
dendrimers. The reciprocal of luq is taken as the time 
unit. The only remaining variable is the free energy, E. 
Simulations were run for values of average ligand num- 
ber n = {0.9, 3.7, 5.8, 13.9}, resulting in a best fit value of 
E = — 3.7± 0.1 kT. The comparison between the monte- 
carlo simulations and the HPLC data are in Figure [2] 
This one-parameter fit to all the distributions is very sat- 
isfactory. 

We also used dendrimers with an average of 68% of 
the active sites blocked by the conjugation of acetamide 
groups (G5 PAMAM; G5(Ac)<7g>(NH2)<34>)- We rep- 
resented this by first using the model above to add ac- 
etamide groups using the acetamide-acetamide catalytic 
interaction free energy barrier Eaa for E in eqn[T] (though 
without steric hindrance; i.e., 5=1). Then we added lig- 
ands, with steric hindrance occurring only between lig- 
ands due to the small size of acetamide. Now there is 
a new parameter, Ea, the ligand-acetamide interaction. 
Thus, the rate for ligand attachment at a site with n; 
ligand neighbors, ria acetamide neighbors, and rif free 
neighbors, where ni + Ua + Uf — 6, is: 



R = LJoL{t)Sini) 



nie 



-E/kT 



+ Uae^^'^/^'^ + nf 



Uf + ni+na 



(2) 



Because the catalyzing amide bond is present in both 
the acetamide- and the ligand-nanoparticle bond, we set 
E = Ea = Eaa- Wc havc no new parameters, but still fit 
the data quite well with one parameter; see Figure [3] 
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Having undersood the distribution of ligands on the 
nanoparticle we now turn to their adsorption and desorp- 
tion from a protein-modified surface. Surface Plasmon 
Resonance (SPR) can sensitively detect the amount of 
material adsorbed onto the surface. In [3] , SPR was used 
to determine the amount of folic acid (FA)-ligated den- 
drimers adsorbed on a surface covered with folate binding 
protein (FBP) as a function of time. This gives the rate 
constants for adsorption, ka, and desorption, k^. It was 
found that k^ for the nanoparticles appeared to decrease 
rapidly with the mean number n of FAs. This was at- 
tributed to multivalent binding. 

Multivalency, i.e. how multiple bonds between ligands 
on the nanoparticle and receptors on the surface effect 
binding is, in general, very complex. However, if C„ is 
broad, apparent multivalent behavior may, in fact, be due 
simply to fluctuations in ligand number, as we will show. 
We explain the data of [3] by supposing that nanoparti- 
cles are permanently bound to the surface if they have 
two or more FA-FBP bonds. Since nanoparticles with no 
FA are not bound at all, only those nanoparticles with 
precisely one FA-FBP bond contribute to kd- Thus for 
the time-scale of the SPR experiments the decrease in 
apparent k^ with increasing n is due to the change in 
ligand distribution C„ and the dilution of the dissociat- 
ing material by the permanently bound material, not an 
enhancement in binding strength. 

This changes the analysis of the SPR data from that 
given in [5]. In the standard multivalent model with 
nanoparticles with different numbers of bonds to the sur- 
face, we expect the dissociation to involve many different 
rates. However, if multiply bound dendrimers do not des- 
orb over the course of the experiement, then there is one 
rate, that of a nanoparticle bound with a single FA-FBP 
bond. We reanalyzed the SPR data (see Figure [4| to get 
new values of this kd- The results (see Figure [5pshow a 
single rate kd ^ 'i ■ 10~^ s~^. 

The number of singly bound nanoparticles is the frac- 
tion of bindable nanoparticles (i.e. with n > 1) that have 
exactly one ligand. Recall that if C„ is Poisson then 
Ci{n) = ne~". However, since unligated nanoparticles 
are not bound at all, the fraction of bindable nanoparti- 
cles with one FA is Fxin) = ne~"/(l — e~"). This simple 
model is adequate at small average ligand number to ex- 
plain the data (see Figure [6]). 

As above, we use kinetic Monte-Carlo to find C„ for FA 
on the dendrimers. From this, we estimate the fraction 
of singly-ligated nanoparticles, and thus what fraction 
of the material we expect to remain bound on the SPR 
surface. The results are shown in figure [6j The SPR 
data is consistent with an approximately constant frac- 
tion of persistent material, whereas our simple model of 
singly-ligated dendrimers declines. This is because the 
number of singly-ligated nanoparticles are an underes- 
timate of the number of singly-bound nanoparticles. A 
nanoparticle might have several ligands, but have only 
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FIG. 4: SPR measurement of nanoparticle coverage. During 
adsorption (rising part of curve), buffer with FA-bound den- 
drimer flows across the SPR chip. During desorption, clean 
buffer flows; kd is obtained by fltting the desorption in the 
second phase, dashed Hue. 
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FIG. 5: Dissociation rate kd versus n, from exponential flts 
to the SPR data with a free horizontal asymptote. The data 
are consistent with a rate independent of fi, fed ~ 3 - 10~^ s~^, 
i.e., the rate of singly-bound nanoparticles dissociating from 
the surface. 



one accessible for binding due to steric effects. Thus the 
theoretical limit is a lower bound for the dissociation, as 
is seen in figure [6j Even without these corrections, the 
simple theory fits the data if n is not too large. 

There are theories of multivalent interactions in 
monodisperse systems in the literature [TDl [TT] . As we 
have pointed out, this sort of treatment is not neces- 
sary if time-scale separation exists between single-ligand 
interactions and multiple-ligand interactions. This is 
the case when thermally released ligand-receptor bonds 
holding the nanoparticle to the surface are likely to re- 
form before the remaining bonds break. Based on the 
folic acid-folate binding protein binding energy estimates 
of Licata and Tkachenko [T^], A£' = \lk^yT , the time 
scale for the desorption of a doubly-bound nanoparticle 
is l/kdeyi-p{/^E/ki,T) « 10^° s « 300 years, a time in- 
accessible to the experiment and to biological processes. 
Notably, the standard proprietary software used to ana- 
lyze SPR data typically assumes that all adsorbed ma- 
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FIG. 6; Fraction of singly-ligated nanoparticles vs. n. The 
points are the fraction of residual material in the SPR exper- 
iment. The crosses are the fraction of singly-ligated nanopar- 
ticles as measured by HPLC. The line is the fraction of singly- 
ligated nanoparticles Fi simulated with the kinetic model 
with parameters E = — 5.0kT, which matched the best fit 
to the HPLC data. The dashed line is _Fi for the Poisson 
case. 



time as the distribution evolves kinetically. 

The model suggests that the strong steric hindrance 
prohibits newly attached ligands having more than one 
neighbor. We have examined our simulation results, and 
we find that this leads to a much larger than chance oc- 
currence of isolated pairs of ligands and, at larger h, lin- 
ear arrangements of ligands. If true, this could have sig- 
nificant biological implications regarding that binding of 
a nanoparticle to a cell, since a ligand is likely to have a 
nearby neighbor. The strong steric hindrance model also 
predicts that saturation of the nanoparticle of around 50 
ligands. Hence the ligand distribution narrows for larger 
n, at some point becoming narrower than a Poisson dis- 
tribution. 

This project has been funded in part with Federal 
funds from the National Cancer Institute, National In- 
stitutes of Health, under Award 1 ROl CAl 19409, and 
with Federal funds from the National Science Foundation 
under Award DMS 0554587. 



terial will eventually desorb. This assumption will result 
in an incorrect measurement in scenarios such as this in 
where there is a separation of time scales. 

Monodisperse systems for which multivalent theories 
are posited may be rarer than anticipated; this is partic- 
ularly the case for those systems in which dispersion is 
meant to be controlled by restricting the free ligand con- 
centration (such as in [71[T31[I1]). Many such multivalent 
systems are in fact composed of polydisperse particles 
with a broad ligand distribution. Using a model such 
as ours, one may deconvolute data obtained from inter- 
esting polydisperse systems, even obtaining estimates for 
quantities such as the free energy of cooperativity for 
ligation. 

Our goal in this paper is twofold: first we present an 
interesting mesoscopic system with its microscopic and 
macroscopic characteristics completely described by a 
simple statistical model. Also, we demonstrate that mul- 
tivalent binding behavior observed in these and related 
chemical systems need not be explained by exotic inter- 
actions, but rather with simple physics combined with 
underlying distribution statistics. Our data gives dis- 
tributions with HPLC and desorption rates from SPR, 
represent a rare case in which high quality information 
is available for both the small and the large scale of a 
mesoscopic system. Our model describes both limits, es- 
timates the free energy of cooperativity of about — 4kT 
for the conjugation of ligands, and can predict ligand 
number distributions for other values of n. This also es- 
tablishes the model as a tool for designing chemical syn- 
theses to attempt to tune the dispersion, for example by 
correctly limiting the initial free ligand concentration or 
by quenching the reaction in progress at the appropriate 
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